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We investigate the spectral features of the emission from a superluminal polarization current whose 
distribution pattern rotates with an angular frequency u) and oscillates with an incommensurate 
frequency fi > u. This type of polarization current is found in recent practical machines designed to 
investigate superluminal emission. Although all of the processes involved are linear, we find that the 
broadband emission contains frequencies that are higher than f2 by a factor of the order of {Q,/lu) 2 . 
This generation of frequencies not required for the creation of the source stems from mathematically 
rigorous consequences of the familiar classical expression for the retarded potential. The results 
suggest practical applications for superluminal polarization currents as broad-band radiofrequency 
and infrared sources. 

PACS numbers: 

I. INTRODUCTION 

Considerable recent interest has been generated by the design g, g, Q] , construction Q and testing Q of novel light 
sources which employ extended distributions of polarization currents moving faster than c, the speed of light in vacuo. 
Although Special Relativity does not allow a charged particle with finite inertial mass to move faster than c, there is 
no such restriction on macroscopic polarization currents, because their superluminally moving distribution patterns 
may be created by the coordinated motion of aggregates of subluminally moving particles Q ||, ||, |h], O, |l2[ [is) . 
Moreover, the distinction between massive and massless charges is not in any way reflected in Maxwell's equations, 
so that the propagating distribution patterns of such polarization currents radiate as would any other moving sources 
of the electromagnetic field |7|, [ll], 12, 13]. 



In view of this interest, the present paper explores the frequency spectrum produced by such superluminal sources. 
We base our analysis on the practical implementation of a superluminal source (Appendix A, p], 0]), and consider 
the electromagnetic waves emitted from a polarization current whose distribution pattern rotates with an angular 
frequency lu and oscillates with an incommensurate frequency Q > lu (Table 1). It is found that the broadband signals 
carried by such waves contain frequencies that are by a factor of the order of (il/ui) 2 higher than the frequency fl; 
i.e. the radiation contains frequencies which are not required for the creation or the practical implementation of the 
source. This results from the cooperation of two effects, both of which are mathematically rigorous consequences of 
the familiar classical expression for the retarded potential. 

1. The retarded time is a multi- valued function of the observation time in the superluminal regime, so that the 
interval of retarded time during which a particular set of wave fronts is emitted by a source point can be 
significantly longer than the interval of observation time during which the same set of wave fronts is received at 
the observation point. 

2. The centripetal acceleration enriches the spectral content of a rotating volume source, for which Q/lu is different 
from an integer, by effectively endowing the distribution of its density with space-time discontinuities. 

These two effects make it possible, for instance, to generate a broadband pulse of radiation whose spectrum has a 
peak in the terahertz band by means of a device whose construction and operation only entails oscillations at two 
radio frequencies: at Q,/(2n) — 302.5 MHz and at a multiple mLu (with, e.g., m = 72) of lu/(2it) — 5 MHz (see Table 
1 and Appendix A). 



This paper is organised as follows. Section [I A outlines the novel features of the radiation from oscillating, acceler- 
ated superluminal polarization currents, contrasting this type of emission with those of fast travelling wave antennas, 
and Cerenkov and synchrotron radiations. Although considerable mathematical rigour is required to derive the fre- 
quency spectrum of such superluminal sources (this stems from the necessary three-dimensionality of superluminal 
sources; superluminal sources cannot be point-like 0, O, |l2|, |l3)), some of the essential features can be demonstrated 
in the case of a localized source, a source whose dimensions are appreciably smaller than the radiation wavelength; 
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Symbol 


Definition 




Angular rotation frequency of the distibution pattern of the source 


n 


The angular frequency with which the source oscillates (in addition to moving) 


f 


Frequency of the radiation generated by the source 


n — 2-irf /lu 


The harmonic number associated with the radiation frequency 


m 


The number of cycles of the sinusoidal wave train representing the azimuthal 
dependence of the rotating source distribution [see Eq. (7)] around the 
circumference of a circle centred on, and normal to, the rotation axis 


|fi ± mu)\ 


The two frequencies to which the spectrum of the spherically decaying component 
of the radiation is limited when Q/uj is an integer. (The non-spherically decaying 
component of the radiation is only emitted at these two frequencies, irrespective 
of whether Q/uj is an integer or not.) 



TABLE I: Definition of the various frequencies and numbers used to describe the source and the emitted radiation 



this is carried out in Section [IB as an aid to understanding the more complex analysis that follows. Section II C 
introduces those properties of a rotating volume source that contribute to to the distinctive spectral features of the 
emitted radiation. The algebraic descripti on of the polarization current (its time dependence and orientation) used in 
the subsequent analys is is given in Section [I D ; all of the later results are presented in terms of the parameters defined 
here. Finally, Section II E emphasizes the important contribution that the centripetal acceleration makes to the broad- 
band nature of the emission. The detailed mathematical formulation of the problem of an extended superluminal 

The special case of those volume elements which contribute coherently to the field at 



source is set out in Section III 



the obse rver is discussed in Section [V , highlighting the importance of focal regions in the space of observation points 
(Section IV A) and the inadequacies of conventional far- field approximations in evaluating the electromagnetic fields 
(Section IV B). The radiation field in the plane of the source's orbit is treated in Section |v|, which contains graphical 
representations of the spectral distribution of the emitted power (Figs. 8 and 9); the predicted frequency spectra show 
broadband, high-frequency emission from a source whose implementation entails only two basic frequencies. A similar 
treatment for the radiation field outside the plane of the source's orbit is given in Section VI, which also describes the 
polarization of the emitted radiation (Table 2). A detaile d co mparison with synchrotron and dipole radiations using 
an analog ous m athematical description is given in Section VII , and the efficiency of the radiative process is estimated 
in Section VIII. Conclusions and a summary are given in Section IX. 



II. DISTINCTIVE FEATURES OF THE EMISSION FROM AN EXTENDED, ROTATING, 
OSCILLATING SUPERLUMINAL POLARIZATION CURRENT 

A. Comparison with fast travelling wave antennas, and Cerenkov and synchrotron radiations 

Superluminally moving surface charges are already encountered in fast travelling wave antennas (m) and certain 
types of leaky waveguides |l5| . The wave fronts emanating from such uniformly moving superluminal sources possess 
an envelope on which the vacuum version of the Cerenkov effect can be observed. However, the source considered in 
the present paper has the time dependence of a travelling wave with both an accelerated superluminal motion and an 
oscillating amplitude; neither of these features are present in extant antennas. 

The superluminal sources considered in the present paper have distinctive features of both the sources of Cerenkov 
and synchrotron radiations: not only do they move with a speed that exceeds the propagation speed of the waves they 
generate (like sources of Cerenkov radiation) but also their motion is centripetally accelerated (like that of sources of 
synchrotron radiation). However, in contrast to sources of these familiar types of emission, they are extended (rather 
than point-like) and they fluctuate in their strength (rather than being time- independent in their own rest frames). 

Acceleration leads to the formation of a cusp in the envelope of the wave fronts that emanate from each volume 
element of the superluminal source a curve along which two sheets of the envelope meet tangentially (Fig. 1). (By 
contrast, the conical envelope of wave fronts in the Cerenkov emission has no cusp.) At any given observation time, the 
radiated field entails a set of such envelopes (each associated with the wave fronts emanating from a specific member 
of a corresponding set of source elements) whose cusps pass through the observation point. These caustics arise from 
those volume elements of the source which approach the observer, along the radiation direction, with the speed of 
light and zero acceleration at the retarded time (Fig. 2). The contribution from the filamentary locus of such source 
elements toward the value of the field at the observation point has certain unexpected properties. Its intensity, for 
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FIG. 1: (a) Envelope of the spherical wave fronts emanating from a source point S which moves with a constant angular velocity 
uj on a circle of radius r — 2.5c/u> (f = ruj/c — 2.5). The circles in broken lines designate the orbit of S and the light cylinder 
rp = c/u) (rp = 1). The curves to which the emitted wave fronts are tangent are the cross sections of the two sheets <f>± of the 
envelope with the plane of source's orbit, (b) Three-dimensional view of the light cylinder and the envelope of wave fronts for 
the same source point S. The tube-like surface constituting the envelope is symmetric with respect to the plane of the orbit. 
The cusp along which the two sheets of this envelope meet touches, and is tangential to, the light cylinder at a point on the 
plane of the source's orbit and spirals around the rotation axis out into the radiation zone (based on M). 




FIG. 2: The bifurcation surface (i.e. the locus of source points that approach the observer along the radiation direction with the 
speed of light at the retarded time) associated with the observation point P for a clockwise source motion. The cusp Cb, along 
which the two sheets of the bifurcation surface meet tangentially, touches the light cylinder (r — 1) at point O. This cusp curve 
is the locus of source points which approach the observer not only with the speed of light but also with zero acceleration along 
the radiation direction. For an observation point in the radiation zone, the spiralling surface that issues from P undergoes a 
large number of turns, in which its two sheets intersect one another, before reaching the light cylinder. Note that the bifurcation 
surface issues from the observation point P and resides in the space (r, ip, z) of source points, while the envelope of wave fronts 
issues from a source point and resides in the space (rp,tpp,Zp) of observation points; the similarity between these two surfaces 
reflects the reciprocity properties of the Green's function for the problem (based on j?j]). 

instance, does not diminish with the distance Rp from the source like Rp~ 2 , as in the case of a spherically spreading 
wave, but more slowly: like Rp~ K with < k < 2 (see Rj). 

In the current paper, we shall restrict our analysis to the frequency content of the spherically-decaying emission; 
the spectrum of the non-spherically-decaying component of the emission is treated in a separate work |p9f. As has 
already been mentioned, a rigorous analysis of this problem is of considerable complexity (Sections 
the extended nature of superluminal sources; we therefore first briefly discuss the case of a localized source, a source 
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whose dimensions are appreciably smaller than the radiation wavelength. This simplified case defines the Green's 
function for the more complex problem, and illustrates some properties of the emission which are also valid in the 
more general case. 



B. Properties of the emission from a localized, rotating and oscillating superluminal source 

Consider a localized source (a source whose dimensions are appreciably smaller than the radiation wavelength, i.e., 
essentially a point source) which moves on a circle of radius r with the constant angular velocity uje z , i.e. whose path 
x(i) is given, in terms of the cylindrical polar coordinates (r, tp, z), by 

r = const., z = const., tp = tp + lot, (1) 

where e z is the basis vector associated with z, and tp the initial value of (p. 

The wave fronts that are emitted by this point source in an empty and unbounded space are described by 

|x P -x(i)| = c(tp-t), (2) 

where the coordinates (xp,ip) = (rp,ipp,Zp,tp) mark the space-time of observation points. The distance |xp— x| = R 
between the observation point xp and the point source x is given by 

R = [{z P - zf + r P 2 + r 2 - 2r P rcos(tpp - tp)]i , (3) 

so that the insertion of Eq. (1) in Eq. (2) yields 

t P = t + [{z P — z) 2 + rp 2 + r 2 - 2r P r cos(tp P — tp — wt)] 5 /c. (4) 

For a given source point (r, tp, z) and various positions (rp, tpp, Zp) of the observation point, this dependence of the 
reception time tp on the emission time t can have one of the generic forms shown in Fig. 3. 

When rui > c, there is a one-dimensional set of observation points for which dR/dt = — c and d 2 R/dt 2 — at the 
emission times t c of the waves, i.e. whose members are approached by the superluminally moving source point with 
the speed of light and zero acceleration at the retarded time. These observation points are located on the cusp of the 
envelope of the emitted wave fronts (Fig. 1), where curve (b) of Fig. 3 passes through an inflection point. In their 
vicinity, Eq. (4) reduces to 

t P =t Pc + ±LU 2 (t-t c ) 3 + ■■■, (5) 

where tp c is the value of tp at which the waves emitted at t = t c arrive, and constructively interfere, at the cusp 
curve of the envelope (see Appendix C of JtJ ) . Note that the coefficient of the third-order term in the above Taylor 
expansion happens to be independent of the coordinates r and z of the source, a feature that enhances the cooperative 
(coherent) nature of the process to be described below (Section |IV|) . 

Now suppose that, in addition to moving faster than its own waves, the point source in question has a strength 
which fluctuates with a frequency higher than that of its rotation, like cos(f2i) with > u> (Table 1). The amplitude 
of the field it would generate will then be proportional to the retarded value of this fluctuating factor, i.e. to a function 
of tp that, according to Eq. (5), has the form cos{0[t c + 6^cu~3 (tp — tp c )^}} at points close to the caustics of the 
wave fronts. The period of oscillations of the observed field is given by the time interval in which the argument of 
this function changes from f2t e to fit c + 2n. It will have a value, therefore, that is by the factor |(7ru;/f2) 2 shorter 
than the period 27r/0 of the fluctuations of the source strength. Stated differently, the time interval t ~ t c — 2tt/Q in 
which a set of wave fronts is emitted is by the factor |(7ru;/f2) _2 longer than the time interval tp — tp c during which 
the same set of wave fronts is received. 



C. Extension to superluminal volume sources 



The field of a uniformly rotating point source constitutes only the Green's function for the present emission pro- 
cess Q. For a corresponding extended source to emit waves of a certain frequency, it is necessary in addition that 
the spectrum of temporal fluctuations of its density should contain that frequency. In this respect, the frequency 
enhancing effect associated with a rotating superluminal source differs radically from that which is familiar from 
synchrotron radiation. The space-time distribution of density for the charged particle which acts as the source of 
synchrotron radiation entails the Dirac delta function and so has a spectral decomposition that is independent of 
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FIG. 3: The relationship between the observation time tp and the emission time t for an observation point that lies (a) inside 
or on, (b) on the cusp of, and (c) outside the envelope of the wave fronts or the bifurcation surface shown in Figs. 1 and 2. This 
relationship is given by g(r,tp, z;rp, tpp, Zp) — tp — tpp, an equation which applies to the envelope when the position (r,tp,z) 
of the source point is fixed and to the bifurcation surface when the location (rp, tpp, zp) of the observer is fixed. (Note that, 
by virtue of the linear relation tp — tp + cut, the motion of the source may be parametrized either by t or by tp.) The maxima 
and minima of curve (a), at which dR/dt = — c, occur on the sheets </>+ an d 4>~ of the envelope (or the bifurcation surface), 
respectively (see Figs. 1 and 2). The inflection points of curve (b), at which d 2 R/dt 2 = 0, occur on the cusp curve of the 
envelope (or the bifurcation surface). 



frequency. That the maximal intensity in the spectrum of synchrotron radiation corresponds to a frequency which is 
much higher than the rotatio n freq uency of its source merely reflects a spectral property of the Green's function for 
that emission process (Section VII). Superluminal sources, on the other hand, are necessarily extended H, [ll], |l2|, [l3| ; 



the spectra of their densities do not, in general, contain all frequencies. 

We now come to a vital distinction between stationary and rotating volume sources which leads to a radical difference 
in the spectral content of the associated emission. This is connected with the differing constraints on the ranges of 
values of tp and tp, constraints which in turn dictate the form assumed by the Fourier decomposition of the source 
density with respect to time. 

The space-time distribution of the rotating point source described in Eq. (1), whose path may be written as r = ro, 
tp = tp + ujt, z — 0, has the density 

p(r,tp,z,t) = qS(r - r )S(tp - wt - tp)S(z)/r, (6) 

where 5 is the Dirac delta function, q is the volume integral of p, and rg a constant. It is known from the analysis 
of synchrotron radiation that the azimuthal angle tp which appears in Eqs. (1) and (6) is not limited to an interval 
of length 2ir, as in the description of a stationary source distribution, but (by virtue of having the value tp + cot) can 
range over the same interval as the time t, i.e. over (—00,00). Since tp = 2n represents the same point in space as 
tp = at any given (r, z), the source density (6) is periodic, both in tp and in t. However, the azimuthal coordinate tp 
does not discontinuously change back to zero each time a rotation is completed. If the time interval over which the 
source density is described by Eq. (6) exceeds a rotation period, then the angle tp = tp + tot that is traversed by the 
source during this time interval would also exceed 2ir. 

Now consider a localized volume source that rotates about the z-axis with a constant angular frequency to (Table 1). 
The density distribution p(r, tp, z, t) for a source of this kind depends on the azimuthal coordinate tp in only the 
combination tp = tp — ut, i.e. is a function of (r, tp, z, t). If we label each volume element of this source by the value 
tp of its azimuthal coordinate tp at t = 0, the equations describing the trajectories of these elements would each have 
the same form as Eq. (1). It can be seen from the collection of space-time trajectories of the constituent volume 
elements of this source, therefore, that the ranges of values of both tp and t are infinite, as in the case of a rotating 
point source, but values of the Lagrangian coordinate tp are limited to an interval of length 2ir, e.g. to — tt < tp < tt 
(Fig. 4). The coordinate tp cannot range over a wider interval because no volume element of an extended source may 
be labelled by more than one value of a Lagrangian coordinate. Phrased differently, the aggregate of volume elements 
that constitute a rotating source in its entirety can at most occupy an azimuthal interval of length 2ir at any given 
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FIG. 4: Projection of the trajectories (world lines) of the volume elements of a uniformly rotating extended source onto the 
(tp,t) space. The Lagrangian coordinate ip = ip — cot designating the initial (t = 0) position of each source element lies in (— n,n) 
while both ip and t range over (— oo,co). The space-time trajectory (world tube) of the extended source itself in (ip,t) space 
consists of the array of trajectories of its constituent volume elements (the broken lines) encompassed by the lines ip = —it and 

(p — TV. 

A 




FIG. 5: The <f> dependence, cos(nvp) cos[f2(</5 — ip)/uo\, of the source density decribed by Eq. (7) at a fixed (r, ip, z) for m = 1, 
fl/uj = 3.5 and <p = n/5. This tp dependence is, by virtue of the relationship t — (tp — ip)/u}, equivalent to a time dependence: the 
space-time of source elements may be marked either by the coordinates (r, tp, z, t), as in Eq. (7), or by the coordinates (r, ip, z, (p) 
as in Eq. (19). The plotted function is physically meaningful only within the interval — -n < ip < n. However, once this function 
is expanded into a Fourier series over the interval (— n.n), a new periodic function results, represented by the series, which is 
coincident with the original function within (— tt,tt) and periodically reproduces the original function outside this interval. The 
periodic function represented by the series outside (— 7r,7r) is designated by broken lines. That the Fourier expansion of source 
density (7) should contain an infinite number of terms may be understood in two (mathematically equivalent) ways: either as 
due to the discontinuities of the periodic function shown here or as due to the non-zero difference between the values of this 
density at p> = ±7r. 

time (e.g. t = 0). 

D. Algebraic representation of polarization currents with superluminally rotating distribution patterns 

For an extended source of radiation whose distribution pattern rotates uniformly, the cylindrical components of 
the electric current density, j, similarly depend on ip only via the Lagrangian coordinate <£> = (p — cut: they have the 
space-time dependence j r ,ip,z(r, <fi, z, t) — j r ,Kp,z{r 1 0, z), where j r ,^p,z stand for the components of j along the cylindrical 
base vectors (e r ,e ¥ ,,e 2 ). In this paper we consider sources that oscillate in addition to moving, i.e. sources for which 
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jr,<p,z are given by jV,^,z(f, <p, z )f(t) with an additional dependence f(t) on time. To be specific, we base the analysis 
on a representative polarization current j = dP/dt for which 

P r , v ,z(r, <p, z,t) = a rtlPtZ (r, z) cos(nvp) cos(fii), —n < ip < n, (7) 

where P r , v ,z & r e the cylindrical components of the polarization (the electric dipole moment per unit volume), s(r, z) 
is an arbitrary vector that vanishes outside a finite region of the (r, z) space and m is a positive integer (Table 1). 

For a fixed value of t, the azimuthal dependence of the above density along each circle of radius r within the source is 
the same as that of a sinusoidal wave train, with the wavelength 2-nr/m, whose m cycles fit around the circumference 
of the circle smoothly. As time elapses, this wave train both propagates around each circle with the velocity rto and 
oscillates in its amplitude with the frequency (Table 1). The vector s is here left arbitrary in order that we may 
later investigate the polarization of the resulting radiation for all possible directions of the emitting current (this 
will be summarised in Table 2). Note that one can construct any distribution with a uniformly rotating pattern, 
Pr,ip,z(r, <p, z), by the superposition over m of terms of the form s rtVtZ (r, z, m) cos(mip). 

An experimentally viable device capable of generating a polarization with the distribution (7) is described in 
Appendix A pL ||. Even though the practical implementation of the polarization (7) by means of this specific device 
only entails the two frequencies mu> and fi, the spectral decomposition of this polarization consists of an infinite set of 
frequencies if Q/lo is different from an integer. One can directly demonstrate this by Fourier analyzing the right-hand 
side of Eq. (7). 

Because the domain of definition of the source density (7) only extends over the interval —tt < ip < n, or equivalently 
(ip — tt)/lj < t < (tp + tt)/w, Fourier decomposition of the time dependence of this function at a given (r, ip, z) should 
be performed by means of a series rather than an integral. Representation of Eq. (7) by a Fourier integral would 
entail assumptions about the dependence of j r ,<p,z on ip in intervals (— oo < ip < —n and 7r < tp < oo) which lie 
outside the domain of definition of this density. The following Fourier series faithfully represents the right-hand side 
of Eq. (7) within its domain of definition and replaces it by a periodic function outside the physically relevant domain 
— 7r < ip < 7r: 

oo 

Pr !ip! z(r,cp,z,t) = s r ,<p, z (r,z) ^2 C n exp(inut), (8) 

n— — oo 

where 

C n = (2tt/uj)^ 1 f(p2n)/u cos{m(ip — wij] cos(fii) exp(— iruvt) 
= (-l) m +"+ 1 (2 7 r)- 1 (n - n/w)[(n - n/cu) 2 - m 2 ]" 1 

x sm(nfl/uj) exp[— i(n — Q/u)(p] + {Q — > — £1}. (9) 

The symbol {tt — > —Q} designates a term like the one preceding it in which f2 is everywhere replaced by —51. 

Note that the Fourier components C n of the function describing the time dependence of P r ,<p,z would vanish for all 
n 7^ (Q/uj) ± m only if VL/uj is an integer (Table 1). When and u> are not commensurable, the factor sin(7rf2/oj) in 
Eq. (9) is different from zero and the spectrum of P r ,ip,z contains all frequencies. It is not difficult to see the reason 
for this: the Fourier series representation of P r ,ip,z equals a periodic function in — oo < tp < oo whose values at the 
beginning and at the end of a period are different when Q/uj is different from an integer (see Fig. 5). The higher 
frequencies stem from these step-like discontinuities of the global function represented by the series, discontinuities 
which lie outside the physically relevant domain — 7r < tp < n but which nevertheless mathematically influence the 
Fourier expansion of the limited part of this function that describes the source density. 

E. The role of centripetal acceleration in providing broadband emission 

Physically, the agent responsible for this implicit discontinuity and so the broadening of the spectrum of the 
source density is centripetal acceleration. That the spectral content of a rectilinearly moving source is broadened by 
acceleration can be directly seen from the transformation x' — x + ut + \at 2 between an accelerated frame S (with 
the initial velocity u and acceleration a) and an inertial frame S' . A source that is spatially monochromatic in its own 
rest frame, e.g. exp(ifcir) with the wavenumber k, is transformed into one, exp[ifc(a/ — ut — ^at 2 )], whose spectrum 
does not even decay at high frequencies. The corresponding effect of centripetal acceleration on the spectrum of a 
source is more subtle and is manifested in a kinematic constraint set by the geometry of rotation. 

Because there is only one parameter (u>) for describing both the speed (rw) and the acceleration (rw 2 ) of a uniformly 
rotating source element, the centripetal acceleration of a rotating extended source shows up in the constraint — n < 
ip < 7r on the range of values of the Lagrangian coordinate (p: an aspect of the geometry of rotation that acts as an 
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additional parameter. Had the range of ip been infinite, the time dependence of the distribution in Eq. (7) would have 
been indistinguishable from that of a rectilinearly moving source with a constant velocity and its Fourier transform 
would have contained only the frequencies |fl ± muo\, irrespective of whether Vl and ui are commensurable or not. It 
is the fact that (p is an angle marking the elements of a rotating source that gives rise to the constraint —n < (p < n, 
to effective space-time discontinuities in the source density, and to the higher frequencies. 

Given that both a representative source density and the Green's function for the present emission process have 
spectra which contain infinite sets of frequencies, it is not unexpected that the analysis in Sections ^ and [Vl| should 
predict a resulting radiation that is correspondingly broadband. This prediction is not incompatible with the fact 
that oscillations at no more than two frequencies (muj and fl) are required for creating the representative source 
distribution [Eq. (7)] in the laboratory (Appendix A; Table 1). As in the case of any other linear system, the present 
emission process generates an output only at those frequencies which are carried both by its input (the source) and its 
response (Green's) function. What entails only two frequencies is the practical implementation of the source we are 
considering and not its spectral content. The density distribution of the present source includes implicit space-time 
discontinuities (Fig. 5) whose Fourier decompositions contain all frequencies [Eq. (9)]. These step-like discontinuities 
do not require (for their practical implementation) the creation of any rapid changes either in space or in time; they 
have to do with the geometry of rotation and automatically stem from centripetal acceleration. 

The radiation field that arises from the superluminal portion (r > c/uj) of the volume source described in Eq. (7) 
consists (as shown in Q) of two components: a spherically decaying component whose intensity diminishes like Rp 
with the distance Rp from the source, and a non-spherically spreading component whose intensity diminishes more 
slowly with distance. The analysis in this paper is concerned only with the spectral properties of the spherically 
decaying component of the radiation. The component of the radiation whose intensity decays like Rp 1 instead of 
Rp~ 2 is emitted only at the two frequencies \Q ± muj\ Q. 



III. DETAILED FORMULATION OF THE PROBLEM 



A. Electromagnetic fields in the far- field limit 



In the absence of boundaries, the retarded potential arising from any localized distribution of charges and 
currents with a density is given by 



A»(x P ,t P ) = c~ x J d 3 xdtf(x,t)S(t P -t-R/c)/R, ^ = 0, • • - , 3 



(10) 



where R stands for the magnitude of R = xp — x, and fi = 1,2,3 designate the spatial components, A and j, of 
and in a Cartesian coordinate system. The expressions that follow for the electromagnetic fields 

E = -VpA° — dA/d(ct P ) and B = V P xA (11) 

when we simply differentiate Eq. (10) under the integral sign, and evaluate the resulting integrals by parts, are 

6{t P - t - R/c) 



E 



d x dt ■ 



R 



ldj 

cdt 



c dt RJ R 



(12) 



B = — - / d 3 xdt 



5{t P -t-R/c) TL fid} j 
R R X \cdt R 



(13) 



since vanishes outside a finite volume. 

Terms of the order of R~ 2 in the above integrands, which do not contribute toward the flux of energy at infinity, 
may be discarded if we are concerned only with the radiation field. Since the problem we will be considering entails 
the formation of caustics, however, we need to treat the phases of the above integrands, i.e. the arguments of the 
delta functions in Eqs. (12) and (13), more accurately. If we replace dj° /d(ct) in Eq. (12) with —V • j from continuity, 
integrate this term by parts and retain only those terms in the integrands of Eqs. (12) and (13) that are of the order 
of R~ x , we obtain 



E 



1 



d 3 xdt 



5(t P -t- R/c) 
R 



AX ( AX |) 



(14) 



and B ~ nxE. Here, we have set the origin of the coordinate system within the source distribution so that |x| <C |xp| 
for an observation point in the far field and R/R can be approximated by the constant vector n = xp/|xp|; the 
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symbol ~ indicates that the expression is valid in the far-field limit. These differ from the standard expressions for 
radiation fields JItJ only in that the argument of the delta function in their integrands is left exact. 

For the purposes of evaluating the integrals in Eq. (14) for the current density that is given by Eq. (7), the space-time 
of source points may be marked either with (x, t) = (r, tp, z, t) or with the coordinates (r, tp, z, t) that naturally appear 
in the description of that rotating source. In fact, once tp is adopted as the coordinate that ranges over (— 7r,7r), the 
retarded position tp of the rotating source point (r, tp, z) as well as the retarded time t could be used as the coordinate 
whose range is unlimited (see Fig. 4). 

The electric current density j = dP/dt that arises from the polarization distribution (7) is given, in terms of tp and 
<P, by 



+{m — > - 



-i(/J,tp - Qtp/uj)](s r e r + s v e v + s z e z ) 

m, — » — fi}, — 7T < tp < 7T, (15) 



where \x± = (fl/uj) ± to and the symbol {m — > 
to and fl are everywhere replaced by — to and 



—m, n — > — f2} designates a term like the one preceding it in which 
-f2, respectively. To put this source density into a form suitable for 



inserting in Eq. (14), we need to express the (^-dependent base vectors (e r ,e,p,e z ) associated with the source point 
(r, tp, z) in terms of the constant base vectors (e rp , e VP , e Zp ) at the observation point [rp, tpp, zp): 











.e z . 





cos(tp — tpp) 
-sw.((p - tp P ) 




sin(<^ — tpp) 
cos(tp — tpp) 
1 













- & Zp _ 



(16) 



Equations (15) and (16) together with the far-field value of n, 



lim n = sin#pe r 

R^oo 



cos t)pe Zf 



= arctan(rp/zp), 



(17) 



yield the following expression for the source term in Eq. (14): 
nx (nxdj/dt) 



exp[-i(fj,tp - fltp/tj)} { [s v cos dp sin(<p - tp P ) 
—s r cosOp cos(tp — tpp) + s z sin 0p]e± + [s r sm(tp — tpp) 

+s v cos{tp — tpp)]e\\ } + {to — > —m, ft — » — f2}, (18) 

where e|| = e VP (which is parallel to the plane of rotation) and ej^ = fixe| comprise a pair of unit vectors normal to 
the radiation direction n. 



B. Green's functions 



Inserting Eq. (18) in Eq. (14) and changing the variables of integration from (x, f) = (r, tp, z, t) to (r, tp, z,ip), we 
obtain 



E 



2 0/ c) 2 Eu=u+ io°° rdr I-^ dz J -J d( P ^ exp(-i/x(p) { (s r G 2 + s^G-Cj 



-[cos9 P (s ip G 2 - s r Gi) + sin 9 P s z G 3 }e^} + {? 



-m,n->-n}, (19) 



where ~ again indicates that the expression is accurate in the far-field limit, and Gi (i = 1,2,3) are the functions 
resulting from the remaining integration with respect to tp: 



G 2 
lG 3 j 



, 8{g - 0) rn i v 

dtp — exp(i\ltpluj) 

R 



C0s((/3 — tpp) 

sin(i^ — tpp) 
1 



(20) 



Here R is as in Eq. (3), <j) stands for tp — tpp with tp = tp — tot and tpp = tpp — u>tp, the function g is defined by 

g = tp - <pp + R, (21) 

with R = Ruj/c, and Atp is the interval of azimuthal angle traversed by the source. The r and z integrations in Eq. 
(19), though extending over the entire (r, z) space, of course receive contributions only from those regions of this space 
in which the source densities s rjVjZ are non-zero. 

Note that it would make no difference to the outcome of the calculation whether one uses the expression in Eq. (19) 
and integrates over the coordinates (tp, tp) as we have done, or one uses Eq. (14) with (x, t) = (r, tp, z, t) and integrates 
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over (t, ip) as is conventionally done. If one follows the conventional procedure and first integrates with respect to 
t, then the constraint —tt < p < tt would show up as the restriction (p — tt)/uj < t < (ip + tt) /to on the range of t 
integration (see Fig. 4). 

The functions Gi(r, p, z; rp, p>p, zp, pp) here act as Green's functions: they describe the fields of uniformly rotating 
point sources with fixed (Lagrangian) coordinates (r, p, z) whose strengths sinusoidally vary with time. The field (19) 
is given by the superposition of the fields of the assembly of such uniformly rotating volume elements from which the 
extended source (15) is built up. In the special case in which SI = 0, i.e. the strength of the source is constant, G3 
reduces to the Green's function called Go in jjj. The singularity structures of Gi are determined by the stationary 
points of the phase function g and so are identical to the singularity structure already outlined in connection with Go 
(see §). 



C. Spectral decomposition of the radiated field 



Spectral decomposition of the radiated field E may be achieved, as in any other time-dependent problem, simply 
by replacing the delta function in Eq. (20) wit h its Fourier representation. Because the integration with respect to p> 
only extends over the interval (— tt, tt) (Section II C ), Fourier decomposition of the p dependence of this delta function 
should be performed by means of a series. 

The integrand in Eq. (19) needs to be faithfully represented only within the range of integration. Representation of 
this integrand by a Fourier integral would entail assumptions about the dependence of j r ,tp,z on <p in intervals which 
lie outside the domain of definition of the source density: in — 00 < p> < —tt and n < p < 00 (see Section |l|). 

Once the delta function S(g — </>) that appears in the integral representation of Gi in Eq. (20) is expanded into a 
Fourier series over the interval — tt < p> < tt, 



5(9 



= (2tt)- 1 J2 exp[-in(g - cf>)}, 



n— — 00 



Eq. (19) becomes 



E = 5r{e + 2^E„ exp(-m^p)}, 

n=l 



(22) 



(23) 



in which 



E,, 



|(w/c) 2 E M=M± Jq°° rdr dz I-n dvv? exp(-i^) {(s r G 2 + JvGi)e|| 



+ [cos9 P (s v G 2 - s r Gi) + sin9 P s z G 3 ]ej_} + {to -> -m,Q -> -ft}, (24) 
and stands for the real part of Z. The functions Gi in this expression are given by 

4(ng - Sl<p/w)], 



Gi = (2tt) 1 exp(in^) / e^/iexp[ 

' Aip 



r/r 




' cos(p> — 


Vp)~ 


h 




sin((y9 — 


<Pp) 


L/s- 




1 





(25) 



with 



(26) 



and constitute the Fourier components of the Green's functions Gi. 

Because the values of (p are limited to an interval of length 27r, the radiation is emitted in harmonics nw of the 
rotation frequency (Table 1). The periodic nature of the motion of the source imposes this constraint despite the fact 
that the source distribution [Eq. (7)] lacks periodicity. [Recall that the ratio fl/uj that appears in the expression for 
the source density, Eq. (7), is different from an integer.] In the regime Q/lo 3> 1, however, the peak of the spectrum 
happens to occur at such a high value of the harmonic number n that this spectrum is essentially continuous. 



IV. LOCI OF COHERENTLY-CONTRIBUTING SOURCE ELEMENTS 
A. The importance of focal regions in the space of observation points 



The filamentary cusps of the envelopes of the wave fronts that emanate from various volume elements of an extended 
superluminal source (Fig. 1) collectively occupy a tubular volume of the xp space, a volume which we shall refer to 
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as the focal region of the space of observation points. At any given observation point P within this focal region, there 
are certain volume elements of the source whose contributions towards the value of the held at the observation time 
tp superpose coherently, i.e. arrive at xp with the same phase. These consist of those elements of the superluminally 
moving source which approach the observer along the radiation direction with the speed of light and zero acceleration 



at the retarded time (Section II B). Or stated mathematically, for large values of the harmonic number n, the 
main contributions towards the value of the multiple integral (24) representing the radiation held E„ come from the 
stationary points of the optical distance tp — t — R/c, given by the function g(r,tp,z) of Eq. (21), that appears in 
the phase of the rapidly oscillating exponential exp(— ing) in the integrand of this integral ]2(| |22j. As a first 
step towards the asymptotic evaluation of the multiple integral in Eq. (24), therefore, we need to identify the loci of 
points at which the derivatives dg/dr, dg/dtp and dg/dz vanish and to expand g into a Taylor series about each of 
its stationary points. 

There is in the present case a point at which all three of the above derivatives are zero. The coordinates of this 
point, which we shall designate as O, are given by 

O : r = 1, tp = tpp + 2n — arccos(l/fp) = tp Ql z — zp, (27) 

where f = rw/c and fp = rpu/c stand for the values of r and rp in units of the light-cylinder radius c/lu. The second 
derivative d 2 g/d(p 2 of g also vanishes at O. Since the derivative of g with respect to tp at fixed (r, tp, z) is proportional 
to the derivative of t — tp + R/c with respect to t at fixed (r,tp,z), the conditions dg/dtp = and d 2 g/ dtp 2 = 
respectively correspond to the conditions dR/dt = — c and d 2 R/dt 2 = which we encountered in Section fll B (see 



Fig. 3). Not only does the point O belong to the locus of source points that approach the observer with the speed of 
light and zero acceleration at the retarded time (i.e. lies on the cusp curve of the observer's bifurcation surface (Figs. 
1 and 2, 0), but it is in fact the point at which this locus touches, and is tangential to, the light cylinder r = c/lu 
(see Fig. 2). 

Whether or not the above stationary point falls within the domain of integration depends on the position of the 
observation point. For a localized source distribution whose dimensions are much smaller than the distance of the 
observer from the source, O would fall within the domain of integration only if the observer lies in the plane of 
rotation, i.e. if the plane which passes through the observation point and is normal to the rotation axis intersects the 
source distribution; otherwise, there would be no source points for which z equals zp. 

For an observer who is located outside the plane of rotation, i.e. whose coordinate zp does not match the coordinate 
z of any source element, only dg/dr and dg/dtp can vanish simultaneously. This occurs along the curve 

C : f = f c (z) = {i(f P + 1) - [\{f 2 P - l) 2 - (z - z P ) 2 ]i}^ 

tp = tpc(z) = tpp + 2ir — arccos(fc/^p). (28) 

In the far-field limit, where the terms (z — zp) 2 /{rp — l) 2 and fc/fp in Eq. (28) are much smaller than unity, this 
curve coincides with the locus 

C b : r=[l + (z-z P ) 2 /(r 2 p -l)]^ 

tp = tpp + 2ir — arccos[l/(ffp)], (29) 

of source points which approach the observer along the radiation direction with the wave speed and zero acceleration 
at the retarded time, i.e. it coincides with the cusp curve of the bifurcation surface (Fig. 2). 



B. The inadequacy of conventional far- field approximations 



The above calculation makes it clear how essential it is that one should start with the exact form of the optical 
distance |xp — x(i)| —c(tp—t) for identifying the loci of its stationary points. The far-field approximation |xp — x(i)| ~ 
|xp|— x-xp/|xp| that is normally introduced at the outset of a calculation in radiation theory |17]] here would obliterate, 
not only significant geometrical features of the loci of these stationary points, but also such determining characteristics 
as the degree of their degeneracy. The far-field approximation would replace the function g by 

g ~ Rp — zcosdp + tp — tpp — fsinOp cos(tp — tpp), Bp » 1, (30) 

where Rp = (rp + zp)i . It is hardly possible to discern the geometrical details of O and C from this expression, let 
alone their nature. 

To preserve the essential features of g about its critical point O, we need to express this function in terms of the 
variables 



£ = tp-cpo, r) = r-l, C = z-z P , 



(31) 
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prior to proceeding to the far-field limit. The resulting exact expression for g(£,?7, C) reduces, when expanded in 
powers of fp 1 , to 

9 = </>o + € - (1 + V) «n £ + [4(1 + ry) sin 2 (£/2) - (1 + t?) 2 sin 2 £ + ry 2 + C 2 ]/(2r P ) + ■ • ■ , (32) 

where </>o = fp + (/?o ~ fp- The term of order fp 1 in this expansion clearly plays a crucial role in determining the 
nature of the stationary point £ = 77 = £ = and so cannot be discarded as in conventional radiation theory For the 
purposes of calculating the asymptotic values of the radiation integrals by the method of stationary phase |2(| |l], ^2| , 
however, it is mathematically permissible to approximate the coefficient of this term by means of a Taylor expansion 
about O. To within the third order in £ and rj, the result is (ry 2 + C 2 )/2 so that 

9 = 4>o + £ - (1 + rj) sine + (ry 2 + C 2 )/(2fp) + • • • . (33) 

This is a more accurate version of the far-field approximation which, in contrast to that appearing in Eq. (30) , exhibits 
the nondegenrate nature of the stationary point O explicitly: neither the coefficient of ry 2 nor that of £ 2 are zero in 
the Taylor expansion of g about O. 

The corresponding expansion of g about a point (fc, <PCi z) on curve C (with an arbitrary coordinate z ^ zp) can 
likewise be found by first rewriting this function in terms of the variables 

p = f-fc, p = p-<pc- (34) 

The result is 

g = <pc - (fp + A* + [&c ~ 2 (! + p/rc)Rc sin^i + 4f c (f c + p) sin 2 (yi/2) + p 2 }^ (35) 
in which we have denoted the value of R on C by 

R c = [{z P - z) 2 + f% - f 2 c ]i (36) 

and have made use of the fact that, according to Eqs. (28) and (36), rc{rp — equals Rc- 

If we now expand the right-hand side of Eq. (35) in powers of R^ 1 (which tends to Rp 1 in the far zone) and 
approximate the coefficient of Rq 1 in the resulting expansion by its value in the vicinity of p = p = 0, we arrive at 

g = 4> c + P - (1 + p/r c ) sin p. + [p 2 + (f 2 c - l)p 2 ]/{2R c ) + • • • , (37) 

where 4>c = Rc + Vc — fp- The term p 2 in the coefficient of R^, 1 in Eq. (37) plays an essential role in specifying the 
degree to which g is stationary on C: p only appears linearly in the earlier terms of the expansion. The term p in 
this coefficient, on the other hand, merely represents a small correction of the order of Rq 1 to the dependence of the 
value of g on p and is of no consequence as far as the asymptotic values of the radiation integrals are concerned. We 
shall therefore neglect the term R^p 2 in Eq. (37) from now on. 

We h ave al ready alluded to the distinction between observation points in and outside the plane of the source 
(Section FVA|) ; we shall now treat these two cases separately in Sections M and [vl. 



V. THE RADIATION FIELD IN THE PLANE OF THE SOURCE'S ORBIT 



A. Treatment of individual volume elements; asymptotic expansion of the Green's function 

Suppose that the observation point is located within the region z < < zp < z > spanned by the orbital planes of 
various volume elements of the source and that its coordinate ipp at the observation time tp is such that the stationary 
point O falls within the range A<p of integration in Eq. (25) (z< < and z> > stand for the extremities of the 
z extent of the source distribution). Then the leading term in the asymptotic expansion of Gi for large Rp and n 
may be obtained by the method of stationary phase: by replacing the phase function g in Eq. (25) with its expanded 
version (33), approximating the coefficient fi of the rapidly oscillating exponential with its limiting value 

.A| =r p - 1 (sinC -1 l) (38) 

at O, and extending the range of integration A(p to (<po — 7r, (fo + 7r) (see [^(J ^l|. 0). 
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The integral that appears in the resulting expression, 

Gi ~ (27T)- 1 exp{in[£ -4> - {v 2 + C 2 )/(2r P )] + iCl Vo /uj} 

x At <%Mo ex P {-i[(n - n/«)^ - nf sin£]}, (39) 

can be evaluated in terms of an Anger function and its derivative: 

Gi rp- 1 exp{in[<£ -fa- (f] 2 + ( 2 )/{2f P )} + Xlxpo/v} 

x[iJ'n-a/u(nr) 3 n -a/u>(nf) - J n -n/u(nf)] , (40) 

where J and J' are Anger's function and the derivative of Anger's function with respect to its argument and the 
symbol ~ denotes asymptotic approximation. The Anger function 3 v (x) is defined by 

J„(X) = (2 7 r)- 1 /:^eexp[-iK-Xsine)] 

— Ju(x) + sin(i/7r) J Q dr exp(—vT — xsinhr), (41) 

in which J v (x) is the Bessel function of the first kind |2^, M . 

There is no difference between an Anger and a Bessel function when v is an integer. The second integral in Eq. 
(41), which constitutes the difference between these two functions when sin(i/7r) ^ 0, has the asymptotic expansion 

poo 

/ drexp(-^r - xsinhr) ~ (1 + v/x)~ 1 X~ 1 - (1 + *Vx)~ 4 X~ 3 

+(9-f/x)(i + iVx)~V 8 + --- C 42 ) 

for large x and positive v (see J24|). The leading contributions toward the values of the Anger functions in Eq. (40), 
therefore, come from the Bessel functions J n -n/ u (nf) and J' n ~n/ui(nf) whose asymptotic values for large n decay 
more slowly than n^ 1 in the superluminal regime. (Here and in what follows, we treat CI as a positive constant.) 

When the argument of J v {x) is smaller than its order and so can be written as x — ^secha for some a > 0, the 
asymptotic values, for large u, of this Bessel function and its derivative are given by 

Jj, (f sech a) ~ (27Ti/tanha)~ 5 exp[i/(tanh a — a)], (43a) 

and 

J' u (v secha) ~ [A-kv/ sinh 2a)~i exp[i/(tanh a — a)] (436) 

(see [^3]). In this regime, the Bessel functions in question decrease exponentially with increasing v. the exponent 
tanha — a is negative for all positive a. But when the argument of J v (x) is greater than its order and can be written 
as x = v sec for some < < ir/2, we have pj| 

J u (i/sec(3) ~ (±m>tan.j3)~% cos[i/(tan/3 - 0) - §], (44a) 

and 

J' v {vsec0) ~ (m// sin 2/3) ~* sm[i/(tan0 - 0) - §]. (446) 

In this case, the Bessel functions in question oscillate with amplitudes that decrease algebraically, like , with 
increasing v. When x and v are equal, these functions decay even more slowly: J v {v) ~ 0.44730 x v~z and J' v (v) ~ 
0.41085 x i/~s (see §|). 

The contributions in Eq. (24) that arise from the source elements in f < 1 — Cl/(nuj), therefore, are exponentially 
smaller than those which arise from f > 1 — 0/ (nuS): the asymptotic values of the Green's functions appearing in Eq. 
(24) are proportional either to J n —u/u( n ^) or to J'n-Q/uiinr) for large n. In particular, the contributions proportional 
to J n -n/uj( n ) and J' n -n/u>( n )y that arise from a source element at the stationary point O, are greater for the terms 
involving Cl in Eq. (24) than for those involving —Cl. 



B. Frequencies of emission from individual volume elements 

The emission is strongest at the frequency with which the source oscillates in addition to moving, i.e. at a value 
of nw for which the integer n is closest to CI/lo. Figures 6 and 7 respectively show the n dependences of the squares 
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FIG. 6: Spectral distribution of the Green's function G2(n) for fi/w = 15.5, normalized by the value G2(16) of this function 
at a harmonic number (n = 16) close to Q/lj. (Frequency / and harmonic number n are related via 2nf — nco.) The inset 
highlights the highest frequency peak of the spectrum. Note that the ranges of frequencies shown in the figure and its inset are 
complementary. This function represents the spectral distribution of the emission arising from the single comoving point O of 
a polarization along the r or the z directions [see Eqs. (24) and (40)]. 
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FIG. 7: Spectral distribution of the Green's function Gi(n), for Q/uj = 15.5, normalized by the value Gi(16) of this function 
at a harmonic number (n = 16) close to fl/u>. (Frequency / and harmonic number n are related via 2nf — ncu.) The inset 
highlights the highest frequency peak of the spectrum. Note that the ranges of frequencies shown in the figure and its inset are 
complementary. This function represents the spectral distribution of the emission arising from the single comoving point O of 
a polarization current which flows in the 95 direction [see Eq. (40) and Eq. (24) for 6p — n/2]. 

of the two Bessel functions J n -n/uj{ n ) an d J'n-n/ui n ) f° r = 15.5, each normalized by its value at n = 16. The 
plotted quantities are respectively proportional to the squares of moduli of the Green's functions Gi and G\ at O [see 
Eq. (40)] . Figure 6 depicts the spectral distribution of the emission that arises from the single comoving point O of 
a polarization along the r or the z directions, and Fig. 7 depicts the spectrum of the corresponding emission from a 
current which flows in the ip direction [see Eq. (24) for 9p = tt/2}. 

The spectral distributions shown in Figs. 6 and 7 confirm the following result that was earlier inferred from time- 
domain considerations in Section ||: the contributions of a rotating point source that is coincident with O, and 
so approaches the observer with the speed of light and zero acceleration at the retarded time, are made over a 
wide range of frequencies and have peaks at harmonic numbers n 3> 1 which are of the order of (fl/uj) 3 when 
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u> <C SI <C nuj. This is a consequence of the fact that in this regime we have J n -n/uj{ n ) ~ (2/n)s Ai[— (2/n)aif2/w] and 
A-fl/uW ~ — (2/n) 3 Ai'[— (2/n) sfi/w] and that the Airy function Ai[— (2/n)3f2/o;] and its derivative peak where 
the magnitude of their arguments is of the order of unity ||23| . 



C. Superposition of the contributions from the constituent volume elements of the source 

The Green's function Gi(r, ip, z; n) calculated in the preceding section represents the contribution to E„ of a specific 
volume element of the uniformly rotating source: that which has the azimuthal coordinate ip at the time t = and 
which moves on a circular orbit of radius r on a plane that is normal to and crosses the rotation axis at z. To find the 
radiation field that arises from an extended source, we must superpose the contributions from the constituent volume 
elements of that source, i.e. we must insert this Green's function in the integral representation (24) of the field and 
perform the integrations with respect to r, ip and z that extend over the localized region of the rest frame occupied 
by that source. 

Equations (24) and (40) jointly yield 

E„ ~ ifp 1 cxp[-i(n0 o - VLp /uj)]Q^ f^dz exp[- finf^ 1 (£ - z P ) 2 ] 

x J °° r df exp[— \mfp 1 (f — 1) 2 ]V + {m — > — m, SI — > -SI}, (45) 

in which 

Q<f> = -|E^ ± dp fi 2 exp[i(n ~ fi)tp] 

= (-l)"+ m sin(7rSl/w)[^ + 2 (n- A i + )- 1 + A i_ 2 (n- M _)- 1 ], (46) 

and V stands for 

V= (s r W 2 + s v W{)e\\ + [cos 6 P (s v W 2 - s r W!) + sm9 PSz W3]e± (47) 

with 

Wi = [iJ'„_n/o)(nf) J n -n/u(nf) -J n -n/w(w)]- (48) 

Here we have replaced r, rp and z by their values in units of the light-cylinder radius (f, fp and z) and have used 
the definition fi± = (SI/lj) ± m to rewrite sin[7r(n — fi>±)] as (— i) n + m + 1 sin(7rfi/cj). The Anger function J and its 
derivative J' are defined in Eq. (41). 

Note that if Sl/o; is an integer, the factors sin[7r(n — /U±)]/(n — fi±) in Eq. (46) would have the value n at n = fi± 
and would vanish at all other n. When SI and ui are not commensurable, on the other hand, neither the numerators 
nor the denominators in these factors can vanish at any n. For non-integral values of n, the quantity Q^, and hence 
the spectrum of the source density (15), is non-zero for all frequencies (see also Section || and Fig. 5). We shall here 
assume not only that the parameter fi/w is different from an integer but also that it is appreciably greater than unity. 
Our interest lies primarily in the high-frequency regime where the radiation frequency nuj is appreciably greater than 
the frequencies \mu> ± Sl| that enter the creation of the source (cf. Appendix A). 

The leading terms in the asymptotic expansions for large n of the integrals over f and z in the expression for E„ can 
both be found by an elementary version of the method of stationary phase [2l| . p2|] . If the lower and upper limits 
of the radial interval in which the source densities s TiH>iZ are non-zero are denoted by f< < 1 and r> > 1, respectively, 
then the asymptotic value of the integral over f in Eq. (45) is given by 

oo 

fdf cxpl-iinrp 1 ^ - 1) 2 ]V - V|. =1 Q r , (49) 
in which 

Q r = j£ dr exv[-\mrp\r - l) 2 ] 

= {wtp/n^iCfo) - C(t ?< ) - i[S( V> ) - %<)]}, (50a) 

with 



V> = [n/(7rr P )]2(f> - 1), 



(506) 
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T7< = [n/(7rrjO] *(*<-!)> ( 50c ) 

and the functions C(r]) and 5 (77) are the Fresnel integrals ^3|. We have integrated over the interval (r<,r>), rather 
than (— 00, 00), in order to obtain expressions that are valid also within the Fresnel distance rp ~ (r> — r<) (nu>)/(7rc) 
from the source. For an observation point that lies at infinity, the integration with respect to f may be directly 
performed over (—00,00), without introducing f< and f>. 

If we now insert Eq. (49) in Eq. (45), and carry out the remaining integration with respect to z in a similar way, 
we arrive at 



E,, 



in which 



with 



and 



^fp exp[-i(n^ - toyo/u)]Q$Q r Q z 'V\ ft=l £ =ip + {m -> -m, Cl -> -O}, (51) 

J. z> dz exp[— iinfp 1 (i — ip) 2 ] 
1 (7rr P /n)i{C(C>) - C(C<) - i[S(C>) - S(C<)]}, (52a) 

C> = [n/(7rfp)]i(5>-5 P ), (526) 
C< = [n/(7rf P )]*(z< -z P ). (52c) 



Here, z> > zp and z< < 2p are the lower and upper limits on the extent of the source distribution in the z direction. 
(The coordinates r and z, without the superscript", are measured in standard units of length rather than in units of 
the light-cylinder radius c/lu.) 

When the z interval z > — z< that is occupied by the source is appreciably smaller than the radius of the light 
cylinder, the observer would have to have a colatitude dp that is quite close to ir/2 for the stationary point O to lie 
within the source distribution. From Eqs. (47) and (48) for Op — n/2, it therefore follows that 

Y \f=i.z=i P = J n-o/ w («)(sroe|| - s z0 e ± ) + iJ' n _n/w(rc)s,poe||, (53) 

in which s r o, v o,zO stand for the values of s riVi2 at the stationary point O. Thus the contributions of the poloidal 
components of the polarization (s r and s z ) to the radiation field are both stronger than, and 90 degrees out of phase 
with, the contribution of the toriodal component s v (see Table 2). 



D. The radiated power and its spectral distribution 

The expression that is given by Eq. (51) for the Fourier component E„ of the radiation field in the plane of rotation 
applies to any frequency nuj for which n 3> 1, provided of course that the distance rp of the observer appreciably 
exceeds both the radius c/lu of the light cylinder and the dimensions r> — r< and z > — z < of the source. If nto is 
also much greater than the frequencies \muj ± f2| that enter the creation of the source, the asymptotic value of the 
radiation field for large frequency reduces to 

E„ - {-lY^fp 1 cxp[-i(n0 o - fl(p /uj)](m 2 + fl 2 /uj 2 ) sm{nn/uj)Q r Q z 

XTi" 1 [J n -n/u(n>)(s r0 e\\ - s z0 e±) + iJ' n -n/u(»)«vO^I|] > ( 54 ) 

for both the additional terms in the Anger functions [Eq. (41)] and the terms associated with — fl [Eq. 51)] would be 
negligibly small. This expression holds not only in the far zone but also within the Fresnel zone. The Fresnel distance 
from the source, which designates the boundary between the near and far zones at a given frequency nto, is defined 
by rp ~ [nLd / (ire)] max{(r> — r<) 2 , (z > — z<) 2 } in the present case. 

When the observation point lies closer to the source than the Fresnel distance, the arguments of the Fresnel integrals 
C and S in Eqs. (50) and (52) are large and so Q r and Q z assume the values 

Q r - (2nf P /n)^ exp(-i7r/4), r P < (n/7r)(f> - f<) 2 , (55a) 



Q z - (27rfp/n)2 exp(-i7r/4), f P < (n/7r)(z> - z<) 2 . (556) 
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n = 2nf/a> 



FIG. 8: Spectral distribution of the radiation that is generated by a poloidal polarization current for fl/w = 15.5 and m = 72. 
The normalization factor is the value Eg7 of the radiation field at a harmonic number (n = 87) close to the peak emission at 
n — (m + Q/lj) or / = (moj + Q)/2n (see Fig. 9); frequency / and harmonic number n are related via 2nf — nu). The plotted 
quantity has the value 1 at n = 87 and decreases monotonically over the range 87 < n < 400 of harmonic numbers not shown 
here (see Fig. 9). The inset highlights the highest frequency peak of the spectrum. 



Beyond the Fresnel distance from the source, on the other hand, the limiting values of C and S are such |23] that 

Q r ~ r>-f<, f P > (n/w)(r> - f<) 2 , (56a) 



Qz~z>-z<, fp > (n/7r)(z> - z<) 2 (566) 

for in the far zone the arguments of the oscillating exponentials in the integrands in Eqs. (50a) and (52a) tend to zero. 
The inclusion of higher order terms in the expansion of g in powers of f — 1 and z — zp would not alter this result: 
the coefficients of these higher order terms contain correspondingly higher powers of fp 1 . 

The radiated power received by an observer in the far zone per harmonic per unit solid angle is given by dP n /dflp = 
cRp 2 |E„| 2 /(27r), where dflp denotes the element sin 6 p d9 p dip p of solid angle in the space of observation points. 
According to Eqs. (54) and (56), this power has the following asymptotic value for high frequency and 9p = ir/2: 

dP n /dn P ~ (2TTc)- 1 ui 2 (m 2 + tt 2 /uj 2 ) 2 sin 2 (7r^/w)(r> - r<) 2 (z> - z<) 2 

X-n- 2 \j„_ n/uJ (n)(s r oen - s z0 e±) +U' n _ n/uJ (n)s v oe\\\ . (57) 

It depends on the harmonic number n like n~ 2 J 2 n ^Qf UJ (n) if the polarization lies in the z or r directions, and like 

n~ 2 J' n-n/uiiri) if the polarization lies in the ip direction. 

Figures 8 and 9 show the emitted power by a poloidal polarization current (s r or s z ) for Q/lo = 15.5 and m = 72. 
The amplitude of this power has in both figures been normalized by its value at the harmonic number closest to the 
source frequency mto + f2, i.e. by the power that is emitted into n — 87. [Note that the peak emission of the device 
described in Appendix A is always close to / = (muj + fl)/2ir; in the case of £l/u> = integer, the only emission occurs 
at this frequency and its companion / = (mw — Q)/2ir. We show elsewhere jlj| that, irrespective of the value of fi/w, 
the non-spherically-decaying part of the emission only contains the frequencies / = |mw ± f2|/27r.] 

The curve shown in Figs. 8 and 9 decreases monotonically from one (the maximum power) at n = 87 to zero at its 
minimum (see Fig. 9). As a result of the presence of the extra factor Q^ 2 ~ n~ 2 , the value of n at which the curve 
in Fig. 8 attains its high-frequency maximum is somewhat lower than that of the corresponding curve shown in the 
inset to Fig. 6. The order of magnitude of this n, however, is still given by (Sl/w) 3 : the result earlier encountered in 
in the context of a point source thus applies also to an extended source. 
Note that, according to Eqs. (54) and (55), the amplitude of E n is independent of the distance rp from the 
source throughout the Fresnel zone: it diminishes like rp -1 only beyond the Fresnel distance, i.e. for rp 3> 
[nu>/(irc)] max {(r> — r<) 2 ,(z> — z<) 2 }. This non-spherical decay of the field amplitude within the Fresnel zone 
is encountered whenever the radiated wave fronts have envelopes; the field of Cerenkov radiation, for instance, decays 



Sections O and V A 
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n = 2%f/a 
350 400 



!faj = 




FIG. 9: Logarithm of the normalized intensity, log 10 |E„/E87| 2 , shown in Fig. 8 over the lower range 30 < n < 400 of harmonic 
numbers (recall that Q/uj — 15.5 and m = 72). Note that the ranges of frequencies in this figure and in Fig. 8 are complementary. 
The contribution from |/i-o->| = muj — Q. (see Table 1) to the radiation field is in this case too small to make a difference to the 
figure even at these lower frequencies. 

like i?p~ 5 within its Fresnel zone. However, the constancy of the amplitude of the radiation with distance is here 
compensated by a steeper dependence of this amplitude on frequency. At the position of an observer who is closer to 
the source than the Fresnel distance, the rate of decay of dP n /d£lp with n is by the factor rT 2 higher than that of 
the power which is radiated in the far zone. 

VI. THE RADIATION FIELD OUTSIDE THE PLANE OF THE SOURCE'S ORBIT 
A. Asymptotic expansion of the Green's function 



When the observation point is located at a colatitude Op which is different from ir/2, the leading contribution to 
the asymptotic value of the radiation held for large Rp and large n comes from those volume elements of the source 
that lie on the curve described in Eq. (28), the curve we designated as C. For Rp ^> 1, C is coincident with the cusp 
curve of the bifurcation surface, Cf,, and the source points in question are those which approach the observer with 
the speed of light and zero acceleration at the retarded time (see Fig. 2). The integrations with respect to cp and r in 
Eqs. (24) and (25) can, as a result, be evaluated by the method of stationary phase once again. 

It can be seen from the far-field limit of Eq. (29) that the cusp curve Cb would intersect the source distribution if 9p 
lies in the interval \0p — ^| < arccos (l/f>), where r> is the radial coordinate of the outer boundary of the source 0. 
Hence, the range of values of Op to which the following analysis is applicable would be as wide as (0, it) if the extent 
r> — c/uj of the supcrlumially moving part of the source is comparable to the radius c/uj of the light cylinder. 

For Op ^ 7r/2, the relevant expansion of the phase function g about the locus of its stationary points is that found 
in Eq. (37). The counterpart of Eq. (39) is therefore given by 



Gi ~ (2n) 1 exp{i[n(^3 - 
x j^d/x/i^exp-t- 



- 0c - ±£ c y 

[(n — Q/uj)/! — 



)+Qtpc/uj}} 
n{f/fc) sin /x]} 



(58) 



in which fc, fc an d Rc, defined in Eqs. (28) and (36), have their far-field values 

rc — CScQp, ipc — <pp + 3tt/2, and Rc — Rp — z cos Op, 

and fi\ c its limiting value 

fA =i? P _1 ( sin n -1 1). 



(59) 
(60) 



[The term Rq 1 fj? in Eq. (37) has, for the reasons given in Section IV, been omitted here.] 
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The integral that appears in Eq. (58) is expressible in terms of an Anger function and its derivative: 

G l Rp' 1 exp{i[n(<£ -\it-R P + zcos9 P - ^Rp 1 p 2 ) + ilip c /w]}Wi (61) 

with 

Wi ee [iJ'n-n/o, (nf sin 9 P ) J n -n/w (nr sin P ) - J„_q/ w {nf sin 9 P )] . (62) 
Hence, Eqs. (24) and (61) jointly yield 

E„ ~ \R P 1 exp{— i[n(R P + §7r) — fl(pc/^}}Q<p dz exp(ini cos#p) 

x J^fdr expl-linRp 1 ^ -esc <9 P ) 2 ]V + {m -m,ft -> -ft}, (63) 

where V differs from the vector V defined in Eq. (47) only in that is everywhere replaced in it by Wi [cf. Eq. 
(62)]. The quantity Q^, resulting from the integration with respect to <p, is the same here as in Eq. (46). 

The dominant contribution towards the asymptotic approximation to the integral over f for high frequency arises 
from the value of its intgrand at the stationary point f = csc9 p of its phase, just as in Eq. (49). The f integral in 
Eq. (63) therefore has the asymptotic value 

/>oo 

jf fdf exp[-im J Rp 1 (f-csc^p) 2 ]V~csc^V|. =csc0p (5 r , (64) 

where 

Q r = (7rfl P /n)i{C(^>) - C(rj<) - i[S(r?>) - S(rj<)]}, (65a) 

with 

*y> = [n/(7r£p)]*(f>-csc0p), (656) 

jy< = [n/(7r£p)]*(f<-csc0p), (65c) 

and f< < csc9 P and f> > csc9 P are the lower and upper limits of the radial interval in which the source densities 
s r ,ip,z are non-zero. 



B. The radiated power and its spectral distribution 

From Eqs. (63) and (64), it now follows that 



E„ ~ \r P exp{-i[n(i? P + |tt) - n<p c /^WQ^QrQ.z + {m -m, ft -> -ft}, (66) 



in which 



with 



= [s r J n -n/w(n) +is v J'n-n/u,(n)]e\\ + [{s v cos9 P 

-s z sm6 P )3 n _Q/ w (n) - is r cos 9 P 3' n _ n/LU (n)] e±, (67) 

/oo 
exp(irii cos 6 p )s TtiPtZ \ c . (68) 
-oo 

The source densities s rjVjZ in s r , v ,z are evaluated along curve C and so are functions of z only. Note that the Anger 
functions that appear in this expression are precisely the same as those appearing in Eq. (53) and depicted in Figs. 6 
and 7: the field E„ depends on 9 P only through the other factors in Eq. (67). 

The corresponding expression for the radiated power per harmonic per unit solid angle is therefore given by 

dP n /dn P ~ (27r)" 1 ccsc 2 0psin 2 (7rft/w)(TO 2 + ft 2 /w 2 ) 2 (r> - r<) 2 n~ 2 |Q z | 2 (69) 
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0P = f 




Sr / 0, S v — S z — 


linear, en, phase= 


elliptic 


S v 0, S r = Sz = 


linear, ey, phase= f 


elliptic 


Sz / 0, S r = S<p = 


linear, ei, phase= ir 


linear, ex 



TABLE II: The state of polarization of the radiation in and out of the plane of rotation for different orientations of the emitting 
polarization current 

when the observation point lies beyond the Fresnel distance Rp ~ (nw/c)(r> — r<) 2 and n appreciably exceeds 
|m±0/w|. Here, as in Eq. (54), both the additional terms in the Anger functions [Eq. (41)] and the terms associated 
with — f2 [Eq. (66)] are negligibly small, so that the Anger functions in |Q Z | may be approximated by Bcsscl functions. 

The rate of decay of |Q Z | with the harmonic number n depends on the length scale of variations of s rjVjZ |c with 
z. The smoother the distribution of these source densities in z, the faster would be the rate of decay of |Q Z | with n. 
In the case of a source with a limited extent in z whose density falls to zero sharply at its z boundaries, such as the 
source described in Appendix A, the quantity |Q Z | decays like rT x for large n. The dependence of the poloidal part 
of the above dP n /dilp on n for |Q Z | ~ nT 1 differs from the corresponding spectrum shown in Fig. 8 only in that its 
amplitude is reduced and its peak is shifted to a slightly lower value of n. 

In the high-frequency regime [where J n -n/u( n ) decreases like n"], the dominant terms of dP n /dQp in Eqs. (57) 
and (69) decay according to a power law n~ a : the index a in this power law has the value 8/3 in the plane of source's 
orbit and a value > 11/3 outside that plane. 

C. Polarization state of the emitted radiation 

Having examined the dependence of the emitted radiation on the polar angle 9p, we finally summarise its po- 
larization and indicate how this is related to the direction s of the emitting current within the device described in 
Appendix A. This summary, for cases in which one of the cylindrical components of s is appreciably larger than the 
other components, is shown in Table 2. 

VII. COMPARISON WITH CERENKOV, SYNCHROTRON AND DIPOLE RADIATIONS 

The features that the above analysis has in common with that of Cerenkov radiation, such as the multi-valuedncss 
of the retarded time and the presence of an envelope of wave fronts at which the phases of the radiation integrals 
are stationary, are easily recognizable. To compare our results with those that are familiar from the analyses of 
synchrotron and dipole radiations, however, we need to give a brief parallel account of certain features of these more 
conventional emission processes in the present notation. 

The electric current density for a uniformly rotating point source (i.e. the source of synchrotron radiation) is 
described by 

j(r, ip, z, t) = quiS(r - r )S(<p - <p )5(z - zo)e v , (70) 

in which the charge q and the coordinates (ro, <po, z ) are all constant. (Recall that (p = ip — uot.) 
Insertion of this source density in Eq. (14) results in 

E ~ q(uj/c) 2 J^'rdr dz J+J dtp6(r- r )6'(tp - <p )5(z - z ) 

x(Gie|| + cos P G 2 e ± ), (71) 

where G\ and G2 are defined by the same integral as that in Eq. (20) with fl = 0. The time dependence of the 
integrand in Eq. (71) can once again be expanded into a Fourier series, not because (p is limited to an interval of 
length 2*7r (here (p has a single fixed value, ip ) but because the source density is periodic. Evaluation of the trivial 
integrals with respect to (r, ip, z) thus results in the following counterpart of Eq. (24): 

E„ ~ -in(w/c) 2 qr (Gie|| + cos 9pG 2 e ± ). (72) 

The definitions of G\ and G 2 in this expression differ from those appearing in Eqs. (25) and (26) only in that in them 
£7 = 0, Aip is (tpp — 7r, ipp + 7r), and (r, ip, z) are everywhere replaced by (ro, (po, zq). 
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Once the phase function g in these definitions is approximated by its far-field value [Eq. (30)] and the integrations 
with respect to tp are performed, we arrive at expressions for Gj whose limiting values for f — > tq are identical to 
the expressions that follow from the far-field versions of Eqs. (61) and (62) when = 0. Insertion of the resulting 
expressions for G\ and G2 in Eq. (72) then leads to 

E„ ~ inq(u}/c) 2 (ro/Rp)exp[in((po - §7r - R P + z cos9 P )] 

x[iJ' n (nfosm0p)ei| + f^ 1 J n (nf sm9p) cot 9pe±)], (73) 

i.e. to the familiar field of synchrotron radiation (cf. |]T^| ). 

The Bessel functions that appear in Eq. (73) have arguments that are smaller than their orders [as in Eq. (43)] 
and so decrease exponentially with increasing n: the speed tquj of the source is subluminal (fo = tqlu/c < 1) in the 
synchrotron process. In contrast, the Bessel functions that appear in the superluminal regime have arguments which 
could equal or exceed their orders and so oscillate with an amplitude that decreases algebraically, like n _ 3 ; , or 
[see Eq. (44) and the paragraph following it]. Were it to exist, a superluminally rotating point source would 
therefore be a much more efficient source of high-frequency radiation than a subluminally rotating one. In fact, as 
can be more directly seen from an analysis in the time domain Q] , the field of a (hypothetical) superluminally moving 
point source is infinitely strong, i.e. has a divergent value, on the envelope of the wave fronts that emanate from it. 

On the other hand, by virtue of being point-like, the source of synchrotron radiation has a spectrum that already 
contains all frequencies. Because the spectral distribution of the source in Eq. (70) is independent of frequency, 
the spectrum of synchrotron radiation is determined solely by the spectral distribution of its Green's function. An 
extended source is radically different in this respect: the spectral content of an extended source of the same type, i.e. a 
rotating source with a density j = j(r, tp, z) whose strength is time independent in its own rest frame, is limited to only 
those wavelengths which characterize the length scales of its variations in tp. Had O been zero for the source described 
in Eq. (7), the spectrum of this source, and hence that of the radiation that arose from it, would have contained only 
the single frequency muj. [What endows the volume source (7) with the broad — albeit rapidly decaying — spectral 
distribution given in Eq. (9) is a new effect, having to do with centripetal acceleration, which would not come into 
play unless the source strength varies with time.] 

In the subluminal regime, volume-distributed charges and currents are typically weak as sources of radiation: the 
contributions from their separate volume elements (those more distant from each other than one radiation wavelength) 
would arrive at the observer with differing phases and so would, as a rule, superpose incoherently. The contributions 
to the radiation field that arise from the source elements located in the vicinity of a stationary point of the optical 
distance g, however, are an exception to this rule. In the case of a superluminally moving extended source, where 
the derivatives of g have zeros, it is possible for the contributions from source elements that are more distant than 
a wavelength to superpose coherently. Not only does the radiation described in Section |y| receive contibutions from 
an interval of retarded time that is by a factor of the order of longer than its period 2n/(nu>), but also the source 
elements that contribute coherently towards this emission occupy r and z intervals that are by the factors ~ and 
~ n longer than the wavelength 27rc/(nw) in the Fresnel and the far zones, respectively [see the values of Q r and Q z 
found in Eqs. (55) and (56)]. 

To compare the present results with those that are familiar from the analysis of dipole radiation, let us consider a 
case in which the polarization P lies in the z direction (i.e. s r = s v = 0) and Vl/u) is an integer. The only non-zero 
values of the quantity appearing in Eq. (46) are in this case those which occur when n = (Q/ui) ± m, so that the 
amplitude of the radiation field E„ reduces to 

|E»| ~ nfp 1 \s\(f > - f<)(z> - z<)n 2 |J„_ 0/w (n)| (74) 

for n = m + fl/u [see Eqs. (45), (54) and (56)]. 

The tp extent of the contributing part of the source is of the order of a wavelength, 2nc/ '(nw). In terms of the total 
dipole moment p = (r> — r < )(z > — z<)[27rc/(no;)]P of the contributing source, therefore, Eq. (74) can be written as 

|E„| ~ irp-^ncj/cfln J n _ n / U (»)p|. (75) 

This differs from the familiar expression for the radiation field of a stationary dipole in the plane normal to its 
direction |ljj by the factor J n -n/u( n )\> a factor which is of the order of ns when n 3> 1. The difference can 
clearly be traced to the phasing of the array of oscillating dipoles that constitute the present moving source (cf. 
Appendix A) : the field calculated in Section [v| receives contibutions from an interval of retarded time that exceeds 
the period 2tt/ (nu>) of the oscillations of its source by the very same factor. 
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VIII. EFFICIENCY OF THE RADIATIVE PROCESS 

The efficiency of the emission process anaiyzed in Sections [v] and VI is essentially independent of the way in which 
a polarization current with a superluminally rotating distribution pattern is created. Our purpose in this section is 
to derive a general expression for estimating the radiation efficiency in the high-frequency regime n ~ (fl/uj) 3 and to 
apply the resulting expression to the particular method of implementing the source density (7) that is described in 
Appendix A. 

If the polarization current density j that acts as the source of the present radiation is produced by the influence 
of an external electric field E cxt on a polarizable medium with electric susceptibility Xe , then the power required for 
maintaining j within a volume V would be Pj n = J v j • E cxt d 3 :E. The induced polarization P is given by x e E oxt , so 
that the polarization current j = dP/dt would have the magnitude |j| ~ x e il|E 0X t|, where is the dominant frequency 
in the spectrum of oscillations of E oxt and hence j. The input power would therefore be of the order of 

P in ~ L r L^L z \i\ 2 /( Xe Q.) (76) 

in terms of |j|, where we have expressed V as the product L r L^L z of the length scales of the source distribution in 
various directions. 

The power that is emitted into the frequency band (An)w centred at nuj ~ {Q,/uj) 2 CI is, according to the analysis 
in Section [v|, given by 

P out ~c- 1 L r 2 L, 2 |j| 2 n- 2 AnAnp ) (77) 

where the solid angle Afip is an estimate of the size of the beam that is emitted into the plane of source's orbit. Note 
that r> — r<, z > — z < and 0|s| in Eq. (57) correspond to L r and L z , and |j|, respectively, and that the high-frequency 
limit of J n (n) has the same order of magnitude as (fl/uj)^ 1 ~ n~3 . 

The above two expressions for the input and output powers imply that the efficiency of the emission process in 
question has a value of the order of 

Pout/Pin ~ Xe{L r /L^){L z uj/c)n-i An AOp, (78) 

in which we have replaced Cl/u by ni . The interval An over which the high-frequency component of the radiation 
is emitted is of the same order of magnitude as n [cf. Figs. 6-9]. But the solid angle Afip in this expression is only 
a small fraction of 4tt (see below): the estimate in Eq. (77) is valid only within the distance L z of the plane of the 
source's orbit. 

Radiation of frequency nuo ~ (CI/lj) 2 CI can be detected also outside the plane of source's orbit. Although the value 
of AJ7p in the corresponding expression for P ou t/Pin is of the order of unity when Op ^ tt/2 (see Section [V]| ), the 
greater steepness of the spectrum reduces the efficiency in this case: the dependence of P ou t /P n on n is by a factor 
of the order of |Q 2 | 2 ~ n~ 2 smaller than that found in Eq. (78) [see Eq. (69)]. 

Consider an experimental device, such as that described in Appendix A, that is built with a polarizable medium 
with the electric susceptibility Xe — 10 an d electrodes that have the dimensions L r ~ ~ 5 cm, L z ~ 1 cm. 
If this device is operated with u/(2tv) ~ 5 MHz and f2/(27r) ~ 300 MHz (the exact value of fl/uj being different 
from an integer), then the efficiency with which the radiation of frequency uuj/(2tt) ~ (£l/oj) 2 £l/(2ir) ~ 1 THz is 
generated by each electrode would be of the order of P ou t/Pin — 3 x 10~ 6 Af2p. Note that (n / V) [L z lo j c) 2 ~ 10 _1 and 
(n/ir)(L r io/c) 2 ~ 1 in this case, so that the boundary between the Fresnel and far zones lies at a distance from the 
source that is shorter than the radius c/u ~ 10 3 cm of the light cylinder [see Eq. (55)]. Beyond the Fresnel zone, the 
radiation that is beamed into the cylindrical region z < < zp < z > surrounding the plane of source's orbit therefore 
subtends a solid angle that is smaller than Afip ~ L z oj/c ~ radians. 



IX. CONCLUDING REMARKS 



We have investigated the spectral features of the intense localized electromagnetic waves that are generated by 
volume polarization currents with superluminally- moving distribution patterns. The analysis is based upon current 
practical devices for investigating emission from accelerated superluminal sources Jl|, [| [| ; these devices (Appendix A) 
produce polarization currents whose distribution patterns rotate and oscillate with two incommensurate frequencies 
(ui and CI). Although the only frequencies entering the production of the emitting currents are raw and f2 (see Table 1), 
we find that the broadband signals from such devices contain frequencies that are higher than the oscillation frequency 
O by a factor of the order of (fi/w) 2 . This does not mean that the linear emission process considered here is capable 
of generating an output at frequencies that are not carried by its input (the source), thus violating the convolution 
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theorem. What is made possible by this process is to generate radiation of a certain frequency from a source whose 
creation does not require that frequency. The spectra of such sources do contain the emitted frequencies [see Eq. (9)]. 

The high frequencies (none of which are required for the practical implementation of the source) stem from the 
cooperation of the following two effects. The retarded time is a multi-valued function of the observation time in the 
superluminal regime, so that the interval of retarded time during which a particular set of wave fronts is emitted by a 
a volume element of the source can be significantly longer than the interval of observation time during which the same 
set of wave fronts is received at the observation point. In addition, a remarkable effect of centripetal acceleration is 
to enrich the spectral content of a rotating volume source, for which il/u> is different from an integer, by effectively 
endowing the distribution of its density with space-time discontinuities. These results are mathematically rigorous 
consequences of the familiar classical expression for the retarded potential. 

The spectral distribution for the emitted radiation is summarised in Figs. 8 and 9, and its possible polarization 
states are listed in Table 2. 

Many features of the radiation that is discussed in the present paper are shared by the high-frequency emissions that 
would arise from superluminally moving sources with generically different trajectories. An example is a superluminal 
source which moves along a straight line with acceleration. In the case of the circularly-moving superluminal source 
considered here, the parameter to determines both the linear velocity of each source element (ruS) and its acceleration 
(ru 2 ). For a superluminal source whose distribution pattern moves rectilinearly, on the other hand, velocity (v) and 
acceleration (a) are two independent parameters. 

The emission (5t) and reception (Stp) time intervals for the waves arising from the source elements that approach 
the observer with the wave speed and zero acceleration are in the rectilinear case, too, related by a cubic equation: 

6t P = \{a/cf[{v/cf - 1]- W + • • • , (79) 

in which v is the retarded value of the source velocity (see Appendix D of [^). Once again, therefore, a moving source 
whose strength fluctuates like cos(ftt) would generate a field which oscillates with a period different from that of its 
source: with the period [27ra/(cft)] 2 [(u/c) 2 - 1]~ 1 (2tt/S1). 

In practice, however, there is a crucial difference between rectilinear and centripetal accelerations. Centripetal 
acceleration, as we have seen, enriches the spectral content of a rotating volume source at the same time as giving 
rise to the formation of caustics and so the compression of dtp relative to St. In contrast, rectilinear acceleration 
requires for its implementation the very frequencies it endows the source with: the range of frequencies with which 
the amplitude of a linearly accelerated source oscillates at a fixed point on its path (within its distribution) is as wide 
as the range of frequencies that feature in the Fourier decomposition of its density with respect to time. 

The rotating superluminal source described in Eq. (7) can be implemented by an experimentally-viable device [|[ |(| 
whose construction and operation only entail oscillations at the two frequencies mu> and f2 (see Appendix A). As a 
source of radiation at frequencies which cannot be normally generated in the laboratory, except by means of large-scale 
facilities such as synchrotrons or free elec tron lasers, the potential practical significance of such a device is clearly 



enormous Jl] ; as we have shown in Section [VTTj, the efficiency is large enough for many spectroscopic applications to 
be viable. 
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XI. APPENDIX A: PRACTICAL IMPLEMENTATION OF THE SOURCE 



The purpose of this appendix is to demonstrate that the polarization described in Eq. (7) can be implemented by 
an experimentally viable device ^] whose construction and operation only entail oscillations at the two frequencies 
mui and Q. 

Consider a circular ring of radius r, made of a dielectric material, with an array of N electrode pairs that are 
placed beside each other around its circumference (Fig. 10). With a sufficiently large value of N (to be determined 
below) , it would be possible to generate a sinusoidal distribution of polarization along the length of the dielectric by 
applying a voltage to each pair independently. The distribution pattern of this polarization can then be animated, 
i.e. set in motion, by energizing the electrodes with time-varying signals. We can synthesize a transverse polarization 
wave cos[m(cp — Lot)] moving around the ring by driving each electrode pair with a sinusoidal signal whose frequency 
is fixed but whose phase depends on the position of the pair around the ring. 
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FIG. 10: View of the experimental device (a) from the side and (b) from above showing an arc of the dielectric medium (dotted 
region), its polarized part (darker dotted region), and the electrode pairs (designated by '±' where on and by '00' where off). 

The frequency mui and the wavelength 2-kt jm of the travelling polarization wave, and hence its speed ru>, can be 
controlled at will by varying the frequency mui with which the electrodes are driven and the phase difference 2wm/N 
between neighboring electrode pairs. Introducing the factor cos(ftt) of Eq. (7) simply corresponds to mixing a second 
frequency into the signal driving the electrodes with a phase that is the same for all electrodes. 

To estimate the required value of N, let us note that the (ip, t) dependence of the polarization that is thus generated 
by the discrete set of electrodes described above has the form 

N-l 

P(cp, t) = cos(fii) n ( fc - Nip/2n) cos[m{ujt - 2nk/N)], (Al) 

k=0 

in which H(x) denotes the rectangle function, a function that is 1 when \x\ < ^ and zero when |x| > \. [For any 
given k, the function H(k — Nip/2ir) is non-zero only over the interval (2k — l)if/N < ip < (2k + l)7r/iV.j When the 
electrodes operate over a time interval exceeding 2tt/uj, the generated current is a periodic function of ip for which 
the range of values of ip correspondingly exceeds the period 27r. 

The Fourier series representation of the function II(fc — Nip/2n) with the period 27r is given by 

oo 

H(k - Nip/2n) = N- 1 + '^2(mr/2y 1 sm(mr/N) cos[n(tp - 2nk/N)}. (A2) 

n=l 

If we now insert Eq. (A2) in Eq. (Al) and use formula (4.3.32) of |2i| to rewrite the product of the two cosines in 
the resulting expression as the sum of two cosines, we obtain two infinite series each involving a single cosine and 
extending over n = 1, 2, • • • , oo. These two infinite series can then be combined (by replacing n in one of them with 
—n everywhere and performing the summation over n = — 1, — 2, • • • , — oo) to arrive at 

P(p, t) = cos(Ot) E^-ooM" 1 sin(mr/A0 

x 12k=o C0S [ muJ t — np + 2n(n — m)k/N], (A3) 

in which the order of summations with respect to n and k have been interchanged and the contribution iV _1 on the 
right-hand side of Eq. (A2) has been incorporated into the n = term: the coefficient (rt7r) -1 sin(n7r/JV) has the value 
N^ 1 when n = 0. 

The finite sum over k can be evaluated by means of the geometric progression. The result, according to formula 
(1.341.3) of |§, is 

iV-l 

cos[mujt — nip + 2n(n — m)k/N] — cos[mujt — nip + ir(n — m)(N — 1)/N] 

k=0 
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x sin[(n — m)ir] csc[(n — m)w/N]. 



The right-hand side of Eq. (A4) vanishes when (n — m)/N is different from an integer. If n = m + IN, where I is an 
integer, the above sum would have the value N cos(mwt — nip), as can be seen by directly inserting n = m + IN in 
the left-hand side of Eq. (A4) . Performing the summation with respect to k in Eq. (A3) , we therefore obtain 



since only those terms of the infinite series survive for which n has the value m + IN with an I that ranges over all 
integers from — oo to oo. 

We have written out the / = term of the series in Eq. (A5) explicitly in order to bring out the following points. 
The parameter N/m, which signifies the number of electrodes within a wavelength of the source distribution, need 
not be large for the factor (mir/N) -1 sm(mn/N) to be close to unity: this factor equals 0.9 even when N/m is 
only 4. Moreover, if the travelling polarization wave cos(m</3) that is associated with the I = term has a phase 
speed rio that is only moderately super luminal, the phase speeds ru>/\l + Nl/m\ of the waves described by all the 
other terms in the series would be subluminal. Not only would these other polarization waves have amplitudes that 
are by the factor |1 + Nl/m]^ 1 smaller than that of the fundamental wave associated with I — 0, but they would 
generate electromagnetic fields whose characteristics (e.g. the peak frequencies of their spectra) are different from 
those generated by the supcrluminally moving polarization wave. 

The fundamental (I = 0) component of the polarization current that is created by the present device, therefore, has 
precisely the same (<p, t) dependence as that which is described in Eq. (7) of the main text. Neither the reduction in 
its amplitude, that arises from the departure of the value of (mir/N) -1 sm(rmr/N) from unity, nor the presence of 
the other lower amplitude waves that are superposed on it make any difference to the fact that the creation of this 
fundamental component only entails the two frequencies muj and fl. Linearity of the emission process ensures that 
the radiation generated by an individual term of the series in Eq. (A5) is not in any way affected by those that are 
generated by the other terms of this series. 

For the distribution pattern of the created polarization current to be moving, it is however essential that the 
number of electrodes per wavelength of this pattern, N/m, should exceed 2. For N/m = 2, the I = — 1 term becomes 
cos[m(ip + cut)] and so represents a wave that has the same amplitude as, and travels with the same speed in the 
opposite direction to, the wave represented by the I = term. In this particular case, the fundamental wave is thus 
turned into a standing wave. 

Note, finally, that the speed of light is easily attainable: if N = 100 electrodes on a circle of radius r = 1 m arc 
driven with the frequency mcu = 600 MHz and the phase difference 2wm/N — 0.125 radians, then the distribution 
pattern of the induced polarization current would consist of m = 2 wavelengths of a sinusoidal wave train which 
would move along the circle with the speed rui = 3 x 10 8 m/s. Moreover, only an arc of such a circularly shaped 
device is needed for generating the present radiation: the pulse that is received at any given observation point arises 
almost exclusively from the limited part of the source which approaches the observer with the speed of light and zero 
acceleration at the retarded time. 
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